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Abstract. 

We consider the quasinormal modes for a class of black hole spacetimes that, 
informally speaking, contain a closely "squeezed" pair of horizons. (This scenario, 
where the relevant observer is presumed to be "trapped" between the horizons, is 
operationally distinct from near-extremal black holes with an external observer.) It 
is shown, by analytical means, that the spacing of the quasinormal frequencies equals 
the surface gravity at the squeezed horizons. Moreover, we can calculate the real part 
of these frequencies provided that the horizons are sufficiently close together (but not 
necessarily degenerate or even "nearly degenerate" ) . The novelty of our analysis (which 
extends a model-specific treatment by Cardoso and Lemos) is that we consider "dirty" 
black holes; that is, the observable portion of the (static and spherically symmetric) 
spacetime is allowed to contain an arbitrary distribution of matter. 
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1. Introduction 

The perturbations of the spacetime outside of a black hole horizon are expected to 
be radiated away, at late times, with a discrete set of complex-valued frequencies; the 
so-called quasinormal mode frequencies of a black hole [1, 2, 3]. Knowledge of these 
modes should have particular importance in gravitational- wave astronomy and, in a 
more speculative scenario, may even provide insight into the very essence of black hole 
entropy (which still lacks a convincing statistical explanation). 

With regard to the latter motivation, an interesting proposal has been put forth by 
Hod [4]. On the basis of Bohr's correspondence principle, Hod has suggested that, in the 
asymptotic limit of a "highly damped" black hole, f the real part of the quasinormal 
frequency should represent a characteristic (transition) frequency for the black hole 
itself. Moreover, it was then argued that the value of this special frequency could be 
used as a means for uniquely fixing the level spacing of the black hole area spectrum. 
(The concept of a uniformly spaced area spectrum for black holes was first advocated 
by Bekenstein [5].) 

To help flesh out these somewhat esoteric statements, let us consider the asymptotic 
behavior of the quasinormal modes for a Schwarzschild black hole. In the case of scalar or 
gravitational perturbations, these asymptotic frequencies are known to take the following 
form [6, 7]: 

+ 0[n- 1 / 2 } as n — > oo , 1 

where n is the quantum number that labels the frequencies and m is the black hole 
mass. § Note that the asymptotic spacing, or the "gap", is just the Schwarzschild 
surface gravity (k s = l/4m). According to Hod [4], one can identify the real part 
of the subleading term, or the "offset", as the "transition frequency" oo so that 
5m = uj = k, s In 3/2-7T . Then, given the Schwarzschild area formula or A = IQnm 2 , it 
immediately follows that 5 A = 327im5m = 4 In 3 equals the spacing between adjacent 
levels in the area spectrum. It is an interesting consequence that, with this particular 
spacing, exp(S') takes on an integer value (where S = A/4 is the black hole entropy) in 
compliance with statistical expectations. 

The above considerations have helped trigger a recent surge of research activity 
in the realm of quasinormal modes. For instance, there has been some discussion on 
linking the intriguing appearance of In 3 with the gauge group of loop quantum gravity 
(e.g. [8, 9]). Furthermore, there has also been a considerable amount of work on 
evaluating the quasinormal frequencies for various black hole models; both numerically 
(e.g. [10, 11]) and analytically (e.g. [12, 13, 14]). For a more thorough list of references, 

| By highly damped, it is meant that the imaginary portion of the mode frequency has become very 
large. This connection follows from the imaginary part being a measure of the inverse relaxation time 
of a radiating black hole. 

§ Here and throughout, all fundamental constants are fixed to unity and a four-dimensional spacetime 
is presumed. 
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see [15] and [16].) 

It is worth pointing out that most of these frequency calculations have been highly 
model specific. (For a recent status report and summary, see [16].) On the other hand, 
it is the opinion of the current authors that significant progress will almost certainly 
require more general considerations. This perspective has prompted a recent study [15] 
in which we investigated the asymptotic mode behavior for a very general class of 
spacetimes: that of "dirty" black holes. These are static, spherically symmetric but 
otherwise generic spacetimes for which there is a central black hole surrounded by 
arbitrary matter fields [17, 18, 19]. (Keep in mind that such geometries also allow for 
multi-horizon scenarios and the presence of a cosmological horizon.) 

Although our prior treatment [15] confirmed that the gap going as the surface 
gravity is a generic result, it was unable to say anything about the real part of the 
offset. In the current paper, however, we are able to rectify this omission for a certain 
class of dirty black hole spacetimes. More specifically, we will now focus on a scenario 
of "squeezed" horizons. That is to say, the observable part of the spacetime is bounded 
by a pair of horizons which are, in a well-defined sense, regarded as being "close". 
(This should not be confused with the case of a nearly extremal black hole where 
the observer is, most typically, regarded as residing in the exterior portion of the 
spacetime.) For this class of spacetimes, it will be shown that the offset can, in fact, 
be calculated exactly in the squeezing limit. This will be accomplished by casting the 
(generic) "scattering potential" into a recognizable, Poschl-Teller form [21]. Notably, 
the same basic technique has been used by Cardoso and Lemos [22] for the specific case 
of a Schwarzschild-de Sitter ("squeezed") geometry; nevertheless, we are now able to 
demonstrate that their results have much more general applicability. 

The remainder of the paper is arranged as follows. In the following section, we 
present the basic formalism and set up conventions. The third section contains the 
main analysis; that is, the calculation of the quasinormal mode frequencies as previously 
discussed. The main body of the paper ends with a brief summary and discussion. 
There is also an appendix where the scattering potential is rigorously derived for axial 
gravitational perturbations in a dirty black hole background (the scalar-perturbation 
calculation is documented in [15]). 

2. The formal preliminaries 

Let us begin by introducing our so-called dirty black hole spacetime. That is, a 
spherically symmetric and static - - but otherwise general - - black hole geometry 
(alternatively, the geometry of a central black hole surrounded by arbitrary matter 
fields). Without loss of generality, the metric for such a spacetime can always be 
expressed as {e.g., [17, 23]) 

-l 

dr 2 + r W, (2) 



ds 1 



_ e -20(r) 



1 - 



2m(r) 





2m(r) 


dt 2 + 
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where <p(r) and m(r) are model-dependent functions (related to the Morris-Thorne 
"redshift" and "shape" functions, respectively [24]). Keep in mind that m(r) retains its 
usual interpretation as a mass parameter, whereas 0(r) measures any deviation (due to 
matter fields) from the "Schwarzschild form" of = —g u ■ 

It often proves useful to have the metric expressed in terms of a generalized tortoise 
coordinate, 



dr* 
dr 



^ _ 2m(r) ' 
r 

so that the line element (2) can now be written as 



ds 2 



-24>(r) 



1 - 



2m{r) 



-dt 2 + dr 2 +r 2 dti 



(3) 



(4) 



Note that our formalism does, in its full generality, allow for spacetimes with more 
than one causal horizon; one of which may be a cosmological horizon (if the spacetime 
is asymptotically de Sitter). For any such horizon — that is, any radius r = which 
satisfies r h = 2m(r- h ) — one can calculate the surface gravity via standard methods [17]. 
More specifically, for a black hole (event) horizon, 



1 d 

2 d7 

1 

e 

2r h 



-4>(r) 



1 - 



<f>( r h) 



1 - 8irp(r h 




(5) 



where p(r) is the energy density. In contrast, one would take the negative of the above 
expression to obtain the surface gravity at a cosmological horizon [25]. In this case, the 
non-negativity of surface gravity immediately implies a lower bound on the radius of any 
cosmological horizon, r c , as a function of the energy density; namely, r 2 > [87rp(r c )] _1 . 
This inequality can be viewed as a generalization of the Nariai bound [26]; that is, the 
upper limit on the mass of a black hole in a Schwarzschild-de Sitter [Kottler] geometry. 

Let us now assume the existence of two non-degenerate horizons; say, r a and r b with 
r& > r a and no other horizons in between. (There can, however, still be other horizons 
in the spacetime, but these would be operationally irrelevant to our treatment — see 
below.) When considering this particular scenario (especially relevant to Section 3), we 
will always view matters from the perspective of an observer who is "trapped" between 
the two horizons. Hence, considerations can safely be restricted to the submanifold 
defined by r a < r < r b . [It is useful to remember that, for this setup, the region 
r G (r a ,rb) maps into the region r* G (— oo, +oo).] Under these circumstances, the 
metric components of interest can always be expressed as 



and 



9r 



\9tt\ 



1 - 



2m(r) (r — r a ) (r b — r) 



h{r) 



-2<Xr) 



1 - 



2m{r) s 



-2<f>(r)'_ 



r a ) (n 



h(r) 



(6) 



(7) 
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where h(r) is some dimensionless function that is well behaved, regular and positive 
throughout the observable portion of the manifold. Note that one can now write the 
mass as 2m(r)/r = 1 — (r — r a )(rb — r)h(r)/r 2 , indicating that h(r) can be physically 
interpreted as (a dimensionless measure of) the energy density of the extraneous or 
"dirty" matter. 

For future convenience, let us re-express some earlier formalism directly in terms 
of the function h(r) or, rather, the modified form (for calculational simplicity) H(r) = 
e~^h(r)/r 2 . The defining expression for the generalized tortoise coordinate (3) now 
becomes 

^ = [(r-r a )(r b -r)H(r)r 1 ; (8) 

whereas the surface gravity (5) can be re-evaluated to give (for the inner and outer 
horizon, respectively) 

K a = ^ (n - r a ) H(r a ) , (9) 

K b = \ (r b - r a ) H(r b ) . (10) 

Since the primary interest of this paper is the calculation of quasinormal mode 
frequencies, let us remind the reader of some pertinent points. One considers small 
perturbations of the spacetime outside of the relevant horizon(s). As is well known 
in the Schwarzschild case, such considerations ultimately lead to a one-dimensional 
Schro dinger- like equation [27]; that is, 

^-V[r(n)}^ = -k 2 ^ , (11) 

where ip = ip [r(r*)] describes the radial behavior of the (massless) perturbation field, 
k is the frequency and, for a scalar (spin j = 0) field in particular, the "scattering 
potential" is found to be 



V(r)=(l-^ 



l(£+l) 2m 

^"2 /y^> 



with i representing the orbital angular momentum {£ — 0, 1, 2, ...). In the case of higher- 
spin (j > 0) fields, there is an often quoted generalization of this result [28] whereby 
one makes the replacement 2m — > 2m(l — j 2 ) in the numerator of the right-most term. 
However, it is certain that this generalization cannot in all circumstances be correct; see 
the appendix for further discussion. 

We have, quite recently [15], generalized the scalar-perturbation potential to the 
generic spacetime described by equations (2) and (4). The one-dimensional Schrodinger- 
like form (11) does indeed persist, although the scattering potential becomes somewhat 
more complicated. More precisely, 

V(r) = e~ 2 ^ ( 1 - (13) 



x 



+ !)_!_ 2 m ( r ) 2 m '( r ) 



^"2 /y^> /y 2 
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where a prime denotes a derivative with respect to r. It is unclear if there is a simple 
way to generalize this result to arbitrary values of spin, j. For instance, the form of the 
potential for the j = 2 case of axial gravitational perturbations (which we derive in the 
appendix) could not have easily been anticipated. (That is, one could not just make 
the type of naive replacement that was discussed above.) Nonetheless, we propose that 
the precise dependence of the potential on j is, qualitatively speaking, inconsequential 
to the prior [15] or the current treatment. 

The quasinormal modes can be regarded as the complex frequency solutions of 
equation (11) when "radiation boundary conditions" are imposed at the boundaries 
of the observable spacetime [1, 3]. (Complex values of k are, in fact, necessitated by 
having a "mostly positive" potential.) With regard to the case of a Schwarzschild black 
hole, the quasinormal mode frequencies are known to be labeled by a discrete quantum 
number [2] and have a well-known asymptotic form, exhibited in equation (1) and amply 
confirmed both numerically [6, 7] and analytically [12, 13]. 

More generally, we have recently shown that, up to the leading asymptotic order, 
the Schwarzschild form for the mode frequencies is, indeed, a generic feature of black 
hole spacetimes [15]. More to the point, the imaginary spacing between the modes (or 
the "gap") is, in the asymptotic limit, exactly equal to the relevant surface gravity under 
quite generic circumstances. Also in the cited work, we have extended considerations 
to an observer who is trapped between two horizons (precisely, the scenario discussed 
earlier in this section) and argued for the following generalization: || 

k qnm {n) = i[n a K a + n b K b ] + 0[1] when n a , n b -> oo . (14) 

Let us re-emphasize, however, that our prior approach [15] - which utilized the first 
Born approximation to the scattering amplitude — could not say anything about the 
real part of the order-unity term or the "offset" . 

3. The main analysis 

We will now proceed to show that, for a certain class of observers, an exact calculation 
of the offset is still "generically" possible. ^[ To elaborate, we will return to our scenario 
of an observer trapped between horizons at r a and r b , and then consider the "squeezed" 
horizon limit of r a — > r b (or vice versa). [This is, of course, just the degenerate or 
extremal horizon limit; nonetheless, we have avoided using such terminology, as the 
usual inference is an observer exterior to both horizons (e.g., an observer outside of 
a Reissner-Nordstrom black hole).] Note that the following analysis generalizes a 
prior work, specific to the case of a Schwarzschild-de Sitter spacetime, by Cardoso 
and Lemos [22]. 

Let us begin here by establishing what is exactly meant by a spacetime with 
squeezed horizons. For the duration, the horizons will be regarded as sufficiently close 

|| The same expression has also been deduced, by independent means, for the specific model of a 
Schwarzschild-de Sitter black hole [29]. However, for a contrary opinion, see [20]. 
By generic, we mean that the spacetime can be arbitrarily dirty. 
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so that A = (r b — r a ) jr a <C 1 . Given the extent of the relevant manifold (r a < r < r b ), 
it immediately follows that, up to corrections of the order A , r a ~ r ~ r b . 

In view of the above, it is clear that much of the prior formalism will simplify. For 
instance, let us recall equation (8) for the tortoise coordinate. Up to corrections of the 
relative order A (always signified by ~), we can write 

^-[(r-rj^-r^r 1 . (15) 

To justify this, recall that H(r) = e-^ r) h(r)/r 2 is a regular, well-defined and positive 
function throughout the observable region [cf, equations (6) and (7)]. We can then safely 
approximate 

H(r a ) ~ H(r) ~ H(r b ) (16) 

in the region of interest r G (r a ,r b ). In particular 

K a ~ K b (17) 

up to corrections of higher order in A. The approximate relation for the tortoise 
coordinate can now readily be integrated to give 

1 



In 



n 



-Lin 



n 



is) 



(r b - r a ) H(r a ) 

which can then be inverted to yield 

r a + r 6 exp [2/c a r*] r a + r b r b - r a 

r — 7z r = — 7: — + — 7, — tanh K a r, . 19 

1 + exp [2ft r*j 2 2 

Applying equation (9), we can directly substitute the above relation for r = r(r*) into 
(r — r a )(r b — r)H(r a ) and obtain (after some manipulations) 

(r - r a ) (r b - r) H{r a ) ~ " (20) 

2 cosh [K a r^\ 

But, given that H{r) ~ H(r a ) throughout the relevant manifold, this also means that 
-Mr) (i 2m(r)\ K a {r b -r a ) 

It is useful to note that the left-hand side is essentially the time-time component of the 
metric [cf, equation (7)], up to a "distortion" which can be regarded as a constant factor 
in the regime of interest. 

The above result is, basically, all that is needed to ascertain the quasinormal 
modes. To demonstrate this, we will initially concentrate on the simplest case of a 
scalar perturbation and then comment on general values of j below. To begin here, 
let us first call upon the generic form of the scattering potential for scalar fields (13). 
Conveniently, one can always re-express this potential in the following compact form [15]: 
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Evaluating the partial derivatives with the help of equation (8), we find that 

V(r) = (r-r a )(r b -r)H(r) (23) 



x 



„-«r)£(£+l) , (r - r a ) (r b - r) H'(r) (r„ + r b - 2r) (r) 



So far, everything is exact; however, applying the "squeezed approximation" (Ac 1 ) 
and equation (21) in particular, we then get 

V(rm) „ e -<Kr a ) «An~r a ) 1(1+1) = 1(1 + 1) 

2 cosh 2 [/t a r*] r\ h(r a ) cosh 2 [re a r*] 

To put it another way, 

V(n) ~ ^ , (25) 
cosh [« a r*] 

where Vo^ is a model-dependent quantity (in principle, always calculable) that is 
constant for any fixed value of t. 

Significantly, equation (25) can be identified as the Poschl-Teller potential [21]; a 
form of potential for which the one-dimensional scattering equation (11) can readily 
be solved. Choosing appropriate plane-wave boundary conditions for the perturbation 
field, \I> ~ exp[ik(t + r*)] as r* — > =foo , one finds that [30] 



1\ Vo,e 1 
n+ 2j + V^-I 



where n = 0,1,2,.... (26) 



This result generalizes equation (20) in reference [22], which was specific to the case of 
a ( "squeezed" ) Schwarzschild-de Sitter black hole. It is easily confirmed that our form 
of the potential (24) reduces to theirs for this special model. 

Take note of the gap or spacing between the levels at large n; this is precisely 
K fl ~ Kb , an outcome which agrees with our prior findings [cf, equation (14)]. Moreover, 
the offset can now be trivially extracted: 




4£±!) _ I (27) 

Besides the surface gravity, the frequencies depend on only one other parameter of the 
spacetime; this being the dimensionless parameter h(r a ). On the basis of dimensional 
arguments, one might expect h(r a ) to be, at the very most, of the order unity and is, 
in fact, exactly unity for the special case of Schwarzschild-de Sitter space. [It is also 
worth pointing out that h(r) depends on the mass parameter, m(r), but not on the 
redshift function, 0(r); that is, the quasinormal modes, even at the level of the offset, 
cannot give us a very detailed account of the dirty matter.] Therefore, one can see that, 
for any £ > 0, the argument of the square root is positive and roughly proportional 
to £ 2 . [Curiously, such explicit £ dependence is notably absent in many studies on the 
quasinormal modes for an exterior observer (e.g., [6, 7, 12, 13]). Nonetheless, there is 
no obvious contradiction here, as no such work caters to our squeezed-horizon scenario.] 
Meanwhile, £ = requires special consideration, since the square-bracket quantity 
in equation (23) for the potential will then vanish up to the order of A. If one goes on 



Dirty black holes: Quasinormal modes for "squeezed" horizons 



9 



to examine the next-order term in the brackets, it becomes evident that the potential 
can no longer be cast into the Poschl- Teller form. Nevertheless, it can be argued that 
this extra factor of A (for £ = the potential is of order A 3 instead of A 2 ) causes the 
potential to become negligible. 

It should be emphasized that, although approximations have been used in attaining 
the Poschl-Teller form, equation (26) provides an exact description of the quasinormal 
mode behavior in the squeezing limit. This is because any neglected terms in the 
potential are of, at most, the relative order A and A oc K a . Hence, we can more 
precisely write 



n+l) + 



(28) 



2) \ h(r a ) 4 

Moreover, the validity of this result does not depend on the horizons becoming 
degenerate or even "nearly degenerate"; they must only be sufficiently close so that 
the condition r& — r a <C r a is satisfied. 

Before concluding, let us briefly comment on the case of general j. If one considers 
axial gravitational (j = 2) perturbations, then the scattering potential takes the form, 
as derived in the appendix, of equation (50). Closely following the prior methodology, 
we find that, in the squeezed-horizon limit, this potential now becomes 

Vir ,) „ = e -,M * (;> ~ r.) W + 1) - 3] . (29) 

cosh [K a r*\ 2 cosh [K a r*\ r A a 

Hence, all of the above outcomes carry through with the simple replacement Vq^ — > ; 
that is, £{£ +1) — > £(£ + 1) — 2 . Although it is difficult to be precise about 
general values of spin, one might be tempted to suggest a "generalization" of the form 
£{£ + 1) — > £{£ + 1) + F[j] , where F[j] is some simple polynomial of leading order 
j 2 . Nevertheless, even though the j dependence probably deviates from such naive 
expectations, it is still clear that: 

—(%) the qualitative features of the analysis will persist, and, 

—(H) for large enough values of £ the dependence on j will become negligible. 



4. Conclusion 



In summary, we have determined the asymptotic frequencies of the quasinormal modes as 
would be measured by an observer who is trapped between a pair of "squeezed" horizons. 
Although these circumstances are arguably special, our analysis is quite general in 
the sense that the spacetime can contain arbitrary matter fields or "dirt". That is, 
we have allowed for the most general horizon geometries in a static and spherically 
symmetric spacetime; thus extending the Schwarzschild-de Sitter treatment of Cardoso 
and Lemos [22]. 

We have found that the "gap" (or asymptotic level spacing) coincides precisely with 
the surface gravity, Kh (at either horizon); in notable agreement with our other recent 
paper [15]. Moreover, the "offset" (or subleading term) was found to go, roughly, as K h £ 
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for large enough £ and n. Let us re-emphasize that our results are exact in the squeezing 
limit, that is, provided that the spacing between the horizons is small in units of horizon 
radius. 

We should also stress that the scenario considered here is distinct from what is 
usually referred to, in the literature, as near extremality or horizon degeneracy. Studies 
on (nearly) extremal black holes usually presuppose an observer who remains exterior to 
both of the horizons. One should not expect the outcomes of this paper to translate into 
this other type of scenario. Indeed, in our related work [15], we have argued that such an 
exterior observer would likely find the quasinormal modes to collapse, asymptotically, 
to zero momentum. + 

Finally, let us ponder upon the meaning (if any) of an offset term that depends so 
strongly on the orbital angular momentum, £, of the perturbation field. To remind the 
reader, the offset of a (for instance) Schwarzschild black hole is believed to asymptote 
towards the same fixed value [/cln(3)/27r] for any (fixed) choice of £ [6, 7, 12, 13]. There 
is, of course, no contradiction here; any scenario with an exterior observer is decidedly 
different from ours. Nonetheless, this "discrepancy" does raise issues with a popular 
conjecture: the asymptotic value of the offset can be employed to fix the spacing between 
black hole area eigenvalues [4, 8]. If this were truly the case, then it seems rather strange 
that the asymptotic value of the offset can, at times, depend on the angular momentum 
of the particles being emitted. That is to say, one would expect this £ independence to 
persist for all types of black holes or, alternatively, why should only certain classes of 
horizons be subject to quantization? It would seem that the status of using quasinormal 
modes, in this particular context, requires further investigation. 
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Appendix: Wave equation for gravitational perturbations 

One can can always see how a gravitational field is perturbed by simply perturbing the 
relevant metric. Holding the stress tensor fixed and studying the linear perturbations, 
one should then be able to extract the gravitational wave equation and, thus, the 
scattering potential for a graviton. Note, however, that this technique says nothing 
about general j = 2 fields; that is, we will be strictly considering j = 2 perturbations 
due to gravity. [In this regard, let us point out that much of the quasinormal literature 
assumes generalizations (from the scalar case) are straightforward, but (as observed by 
Cvetic and Larsen, and by Kanti and March-Russell [31]) such a naive approach fails 
to incorporate the Bianchi identities except in the vacuum equations. Which is to say, 

+ Actually, the extremal situation is somewhat more complicated than this. See [15] for a detailed 
discussion. 
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the wave equation should, generally, depend on the type of field being discussed and not 
just its spin.] The computation reported below is an independent consistency check on 
that of Karlovini [32] . 

Allowing for the possibility of gravitational radiation, we will (closely following 
Chandrasekhar [33]) consider small perturbations of the metric and characterize these 
by q r ,q t and q e : 

-l 



ds 2 



2m r ) 



dt 2 + 1 



2m r 



r 



dr 2 + rW 



(30) 



+r 2 sin 2 9 [dy? - q r (r, 6, t)dr - q 9 {r, 6, t)d0 - q t {r, 6, t)dt} 2 . 

(So far, all the functional dependence is explicit in the metric but, for the sake of 
brevity, this will not always be so in the following.) The above equation describes the 
so-called axial perturbations, for which a change in the sign of (p necessitates that the 
perturbations must also change sign (for the metric to remain unchanged). It should 
be emphasized that these perturbations are not completely general; one could, just as 
well, always introduce small perturbations into the functions m{r) and <f>(r). This other 
(neglected) type are called polar perturbations; for this kind, a reversal of f does not 
affect the metric. In the case of a vacuum, axial and polar perturbations must decouple 
and any arbitrary perturbation can be written as a linear superposition. Let us note, 
however, that even for the simple model of a Schwarzschild black hole in a vacuum, the 
scattering potential for polar perturbations does not generalize to the (naive) suggestion 
of replacing 2m with 2m(l — j 2 ) in equation (12) for scalars. 

We now proceed by demanding that the perturbations do not change the stress- 
energy tensor (at least) to first order in q. The non-vanishing first-order components of 
the Ricci tensor are found to be 



5R 



13 



r 3 ^Jl - 2m/r sin 2 9 e~^> 



q [rfi \r 2 sin 3 i 



1 - 



t) 



,-<t> 



(31) 



r sin 6 



-20 



yjl - 2m/', 



5R 



23 



(r 2 sin 2 i 



g[ ri 0]r 2 sin 3 6> 

sin# 

e -^(\ - 2m/r 



= Q[t,r],t , 

L 2m^j 

-Q[t,e],t , 



(32) 



5R 



1 



03 



" 3 sin 6 



1 - 2m/r q [t . 



A sin 3 fl N 



(33) 



re~ 



>y/l - 2m/r 



(q[t,e]r 3 sin 3 



which yield the relations 
SR 13 = 



: sin 3 6 1 



2m \ 



r 4 sin 3 6 



-V[t,r],t 



(34) 
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5R 



23 



<=> 



q[ tfi f sin 0(1 —j e~ 



r 2 sin 3 i 



12 

Q[t,e],t , (35) 



e~*(l - 2m/r 

with the third possible expression being redundant. These equations can be most easily 
found by tetrad methods or, alternatively, consult §24 of Chandrasekhar [33] [equations 
(11) and (12)]. Note that these formulas even hold true in the more general case of g u , 
g rr and gee depending on as well. 

To help simplify the above results, let us define the variable Q by 

2m{r) 



Q(r,9,t) = r 2 sm d 9 



(36) 



The equations governing the perturbations [(34 and 35)] can now be expressed in the 
following form: 



r 2 sin 3 

= 7^71 /„\ 9[t,6],t 



(37) 
(38) 



e~*(l - 2m/ r) 
Q, e = -r 4 sin 3 0e^ q [t , r] , t . 

Next, we can solve for separable harmonic perturbations by imposing a suitable 

ansatz, 

Q(r,9,t) =rR(r) tf(0) e ikt , (39) 

with other solutions obtainable by superposition. Also applying the commutativity 
property of the partial derivatives, we can eliminate the function q t from equations (37) 
and (38) to obtain 



2e~*(l - 2m/ r) 



r 2 sin 3 



(rR), 



+ 



Re~ 



sin 3 



(40) 



k 2 {qe,r - q r ,e)e 



—ikt 



2k 2 R'd 



r sin 0(1 — 2m/r)e~^ ' 



which implies 



~R~ 



e _ *(l - 2m/r) 



+r 2 e 2 ^ 



(ri2) >r 
k 2 



+ 



sin 3 (0) 
. 



sin 3 (0) 



(*), 



(41) 



(1 — 2m/r) 

The terms are now in a form that is particularly suitable for finding separable solutions. 
Denoting the separation constant by 7 2 , one finds that 



rV 



sin 3 



- 2m/r) 



rR), 



+ r 2 e 2< ^- 



A: 2 



1 — 2m/r) 



sin 3 



R - 7 2 i? = , 



= . 



The angular equation is quite easily solved and gives 
(0) = sin 2 6 \A P 2 £ (cos(0)) + B Q e 2 (cos(0))~ 



(42) 



(43) 



(44) 
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where 



V9 + 4 7 2 - 1 



(45) 



and Q m are the Legendre polynomials, while A and B are arbitrary constants. To 
ensure that the solutions are real, we must restrict the values of t to integers, and so 

^ = f + £-2 = (£ + 2)(£-l) where £ = 0,1,2,.... (46) 

The radial equation is not so trivially solved. However, we can still transform it 
into a recognizable form by first taking note of equation (3) for the (generalized) tortoise 
coordinate. It follows that the radial equation can also be expressed as 



or 



d 2 R 
dr? 



dr R R ,, 
+ — 



dr 



+ rV 



7 2 r, . 2 



dr 



dr 
dr* 




k 2 R - YR = 




-k 2 R 



(47) 



(48) 



The above can readily be identified as having the form of a one-dimensional 
scattering equation [cf, equation (11)] with a potential of 



or 



dr 



V{r) = e~ 2<t) (l 




2m ^ 



+ 1) 



6m 2m r 



i _ 2m 



(49) 



(50) 



(As a consistency check, it can be verified that this form correctly reduces to the 
potential given in [22] for the same type of perturbations in a Schwarzschild-de 
Sitter background.) This equation can be somewhat simplified by noting that the 
two derivative terms appearing above can be rewritten in terms of the Rg§ (or R^) 
orthonormal component of the Ricci tensor: 



V(r) 



-2<f> 



l - 



2m\ 
r J 



+ 6m(r) 



Equivalently, in terms of the Einstein tensor 

2m\ 



V(r) 



3 -2</> 



r / 



1(1 + 1) _ Mr) G u - a, 



(51) 



(52) 



r° 2 

If desired, the Einstein equations could the be used to rewrite the trailing term in terms 
of the density and radial pressure. Comparing this potential with equation (13) for a 
scalar perturbation, one can see that there is no obvious generalization when the black 
hole is truly "dirty". As a further point of interest, one can also see [most easily from 
equation (49)] that, in the limit of squeezed horizons, only the r j 2 term will contribute 
to the scattering. 

As a second consistency check this can be compared with the analysis of 
Karlovini [32], to which it is equivalent for j = 2. In fact Karlovini also performed 
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the j = I calculation for the Maxwell field. Combined with the (minimally coupled) 
scalar result as discussed earlier in this article it is possible to re-cast Karlovini's result 
in the form 



V(r) = e~ 2 * (l 



2m\ 
r ) 



+ 1) 



+ (i-i 2 ) 



2m{r) 



-(i-j)R§ 



(53) 



for j = 0, 1, 2. Given the results of the scalar, Maxwell, and gravity computations, 
deriving this is a curve-fitting exercise on a curve with 3 points. There is therefore no 
particular reason to trust this formula for higher values of j. 
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